3.7.71 \(\int (a+b \sin (e+f x)) (c+d \sin (e+f x))^3 \, dx\) [671]

Optimal. Leaf size=171 \[ \frac {1}{8} \left (8 a c^3+12 b c^2 d+12 a c d^2+3 b d^3\right ) x-\frac {\left (4 a d \left (4 c^2+d^2\right )+3 b \left (c^3+4 c d^2\right )\right ) \cos (e+f x)}{6 f}-\frac {d \left (6 b c^2+20 a c d+9 b d^2\right ) \cos (e+f x) \sin (e+f x)}{24 f}-\frac {(3 b c+4 a d) \cos (e+f x) (c+d \sin (e+f x))^2}{12 f}-\frac {b \cos (e+f x) (c+d \sin (e+f x))^3}{4 f} \]

[Out]

1/8*(8*a*c^3+12*a*c*d^2+12*b*c^2*d+3*b*d^3)*x-1/6*(4*a*d*(4*c^2+d^2)+3*b*(c^3+4*c*d^2))*cos(f*x+e)/f-1/24*d*(2
0*a*c*d+6*b*c^2+9*b*d^2)*cos(f*x+e)*sin(f*x+e)/f-1/12*(4*a*d+3*b*c)*cos(f*x+e)*(c+d*sin(f*x+e))^2/f-1/4*b*cos(
f*x+e)*(c+d*sin(f*x+e))^3/f

________________________________________________________________________________________

Rubi [A]
time = 0.16, antiderivative size = 171, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, integrand size = 23, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.087, Rules used = {2832, 2813} \begin {gather*} -\frac {d \left (20 a c d+6 b c^2+9 b d^2\right ) \sin (e+f x) \cos (e+f x)}{24 f}-\frac {\left (4 a d \left (4 c^2+d^2\right )+3 b \left (c^3+4 c d^2\right )\right ) \cos (e+f x)}{6 f}+\frac {1}{8} x \left (8 a c^3+12 a c d^2+12 b c^2 d+3 b d^3\right )-\frac {(4 a d+3 b c) \cos (e+f x) (c+d \sin (e+f x))^2}{12 f}-\frac {b \cos (e+f x) (c+d \sin (e+f x))^3}{4 f} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[(a + b*Sin[e + f*x])*(c + d*Sin[e + f*x])^3,x]

[Out]

((8*a*c^3 + 12*b*c^2*d + 12*a*c*d^2 + 3*b*d^3)*x)/8 - ((4*a*d*(4*c^2 + d^2) + 3*b*(c^3 + 4*c*d^2))*Cos[e + f*x
])/(6*f) - (d*(6*b*c^2 + 20*a*c*d + 9*b*d^2)*Cos[e + f*x]*Sin[e + f*x])/(24*f) - ((3*b*c + 4*a*d)*Cos[e + f*x]
*(c + d*Sin[e + f*x])^2)/(12*f) - (b*Cos[e + f*x]*(c + d*Sin[e + f*x])^3)/(4*f)

Rule 2813

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[(2*a*c +
b*d)*(x/2), x] + (-Simp[(b*c + a*d)*(Cos[e + f*x]/f), x] - Simp[b*d*Cos[e + f*x]*(Sin[e + f*x]/(2*f)), x]) /;
FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0]

Rule 2832

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[(-d
)*Cos[e + f*x]*((a + b*Sin[e + f*x])^m/(f*(m + 1))), x] + Dist[1/(m + 1), Int[(a + b*Sin[e + f*x])^(m - 1)*Sim
p[b*d*m + a*c*(m + 1) + (a*d*m + b*c*(m + 1))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[
b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && GtQ[m, 0] && IntegerQ[2*m]

Rubi steps

\begin {align*} \int (a+b \sin (e+f x)) (c+d \sin (e+f x))^3 \, dx &=-\frac {b \cos (e+f x) (c+d \sin (e+f x))^3}{4 f}+\frac {1}{4} \int (c+d \sin (e+f x))^2 (4 a c+3 b d+(3 b c+4 a d) \sin (e+f x)) \, dx\\ &=-\frac {(3 b c+4 a d) \cos (e+f x) (c+d \sin (e+f x))^2}{12 f}-\frac {b \cos (e+f x) (c+d \sin (e+f x))^3}{4 f}+\frac {1}{12} \int (c+d \sin (e+f x)) \left (12 a c^2+15 b c d+8 a d^2+\left (6 b c^2+20 a c d+9 b d^2\right ) \sin (e+f x)\right ) \, dx\\ &=\frac {1}{8} \left (8 a c^3+12 b c^2 d+12 a c d^2+3 b d^3\right ) x-\frac {\left (4 a d \left (4 c^2+d^2\right )+3 b \left (c^3+4 c d^2\right )\right ) \cos (e+f x)}{6 f}-\frac {d \left (6 b c^2+20 a c d+9 b d^2\right ) \cos (e+f x) \sin (e+f x)}{24 f}-\frac {(3 b c+4 a d) \cos (e+f x) (c+d \sin (e+f x))^2}{12 f}-\frac {b \cos (e+f x) (c+d \sin (e+f x))^3}{4 f}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 0.80, size = 143, normalized size = 0.84 \begin {gather*} \frac {-24 \left (3 a d \left (4 c^2+d^2\right )+b \left (4 c^3+9 c d^2\right )\right ) \cos (e+f x)+8 d^2 (3 b c+a d) \cos (3 (e+f x))+3 \left (4 \left (8 a c^3+12 b c^2 d+12 a c d^2+3 b d^3\right ) (e+f x)-8 d \left (3 a c d+b \left (3 c^2+d^2\right )\right ) \sin (2 (e+f x))+b d^3 \sin (4 (e+f x))\right )}{96 f} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[(a + b*Sin[e + f*x])*(c + d*Sin[e + f*x])^3,x]

[Out]

(-24*(3*a*d*(4*c^2 + d^2) + b*(4*c^3 + 9*c*d^2))*Cos[e + f*x] + 8*d^2*(3*b*c + a*d)*Cos[3*(e + f*x)] + 3*(4*(8
*a*c^3 + 12*b*c^2*d + 12*a*c*d^2 + 3*b*d^3)*(e + f*x) - 8*d*(3*a*c*d + b*(3*c^2 + d^2))*Sin[2*(e + f*x)] + b*d
^3*Sin[4*(e + f*x)]))/(96*f)

________________________________________________________________________________________

Maple [A]
time = 0.32, size = 182, normalized size = 1.06

method	result	size
derivativedivides	\(\frac {a \,c^{3} \left (f x +e \right )-3 c^{2} d a \cos \left (f x +e \right )+3 a c \,d^{2} \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-\frac {d^{3} a \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )}{3}-c^{3} b \cos \left (f x +e \right )+3 b \,c^{2} d \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-c \,d^{2} b \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )+b \,d^{3} \left (-\frac {\left (\sin ^{3}\left (f x +e \right )+\frac {3 \sin \left (f x +e \right )}{2}\right ) \cos \left (f x +e \right )}{4}+\frac {3 f x}{8}+\frac {3 e}{8}\right )}{f}\)	\(182\)
default	\(\frac {a \,c^{3} \left (f x +e \right )-3 c^{2} d a \cos \left (f x +e \right )+3 a c \,d^{2} \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-\frac {d^{3} a \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )}{3}-c^{3} b \cos \left (f x +e \right )+3 b \,c^{2} d \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-c \,d^{2} b \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )+b \,d^{3} \left (-\frac {\left (\sin ^{3}\left (f x +e \right )+\frac {3 \sin \left (f x +e \right )}{2}\right ) \cos \left (f x +e \right )}{4}+\frac {3 f x}{8}+\frac {3 e}{8}\right )}{f}\)	\(182\)
risch	\(a \,c^{3} x +\frac {3 x a c \,d^{2}}{2}+\frac {3 x b \,c^{2} d}{2}+\frac {3 x b \,d^{3}}{8}-\frac {3 \cos \left (f x +e \right ) c^{2} d a}{f}-\frac {3 \cos \left (f x +e \right ) d^{3} a}{4 f}-\frac {\cos \left (f x +e \right ) c^{3} b}{f}-\frac {9 \cos \left (f x +e \right ) c \,d^{2} b}{4 f}+\frac {b \,d^{3} \sin \left (4 f x +4 e \right )}{32 f}+\frac {d^{3} \cos \left (3 f x +3 e \right ) a}{12 f}+\frac {d^{2} \cos \left (3 f x +3 e \right ) b c}{4 f}-\frac {3 \sin \left (2 f x +2 e \right ) a c \,d^{2}}{4 f}-\frac {3 \sin \left (2 f x +2 e \right ) b \,c^{2} d}{4 f}-\frac {\sin \left (2 f x +2 e \right ) b \,d^{3}}{4 f}\)	\(204\)
norman	\(\frac {\left (a \,c^{3}+\frac {3}{2} a c \,d^{2}+\frac {3}{2} b \,c^{2} d +\frac {3}{8} b \,d^{3}\right ) x +\left (a \,c^{3}+\frac {3}{2} a c \,d^{2}+\frac {3}{2} b \,c^{2} d +\frac {3}{8} b \,d^{3}\right ) x \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\left (4 a \,c^{3}+6 a c \,d^{2}+6 b \,c^{2} d +\frac {3}{2} b \,d^{3}\right ) x \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\left (4 a \,c^{3}+6 a c \,d^{2}+6 b \,c^{2} d +\frac {3}{2} b \,d^{3}\right ) x \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\left (6 a \,c^{3}+9 a c \,d^{2}+9 b \,c^{2} d +\frac {9}{4} b \,d^{3}\right ) x \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )-\frac {18 c^{2} d a +4 d^{3} a +6 c^{3} b +12 c \,d^{2} b}{3 f}-\frac {2 \left (3 c^{2} d a +c^{3} b \right ) \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{f}-\frac {2 \left (9 c^{2} d a +2 d^{3} a +3 c^{3} b +6 c \,d^{2} b \right ) \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{f}-\frac {2 \left (27 c^{2} d a +8 d^{3} a +9 c^{3} b +24 c \,d^{2} b \right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 f}-\frac {3 d \left (4 a c d +4 b \,c^{2}+b \,d^{2}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{4 f}+\frac {3 d \left (4 a c d +4 b \,c^{2}+b \,d^{2}\right ) \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}-\frac {d \left (12 a c d +12 b \,c^{2}+11 b \,d^{2}\right ) \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}+\frac {d \left (12 a c d +12 b \,c^{2}+11 b \,d^{2}\right ) \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{4}}\)	\(488\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+b*sin(f*x+e))*(c+d*sin(f*x+e))^3,x,method=_RETURNVERBOSE)

[Out]

1/f*(a*c^3*(f*x+e)-3*c^2*d*a*cos(f*x+e)+3*a*c*d^2*(-1/2*cos(f*x+e)*sin(f*x+e)+1/2*f*x+1/2*e)-1/3*d^3*a*(2+sin(
f*x+e)^2)*cos(f*x+e)-c^3*b*cos(f*x+e)+3*b*c^2*d*(-1/2*cos(f*x+e)*sin(f*x+e)+1/2*f*x+1/2*e)-c*d^2*b*(2+sin(f*x+
e)^2)*cos(f*x+e)+b*d^3*(-1/4*(sin(f*x+e)^3+3/2*sin(f*x+e))*cos(f*x+e)+3/8*f*x+3/8*e))

________________________________________________________________________________________

Maxima [A]
time = 0.31, size = 189, normalized size = 1.11 \begin {gather*} \frac {96 \, {\left (f x + e\right )} a c^{3} + 72 \, {\left (2 \, f x + 2 \, e - \sin \left (2 \, f x + 2 \, e\right )\right )} b c^{2} d + 72 \, {\left (2 \, f x + 2 \, e - \sin \left (2 \, f x + 2 \, e\right )\right )} a c d^{2} + 96 \, {\left (\cos \left (f x + e\right )^{3} - 3 \, \cos \left (f x + e\right )\right )} b c d^{2} + 32 \, {\left (\cos \left (f x + e\right )^{3} - 3 \, \cos \left (f x + e\right )\right )} a d^{3} + 3 \, {\left (12 \, f x + 12 \, e + \sin \left (4 \, f x + 4 \, e\right ) - 8 \, \sin \left (2 \, f x + 2 \, e\right )\right )} b d^{3} - 96 \, b c^{3} \cos \left (f x + e\right ) - 288 \, a c^{2} d \cos \left (f x + e\right )}{96 \, f} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*sin(f*x+e))*(c+d*sin(f*x+e))^3,x, algorithm="maxima")

[Out]

1/96*(96*(f*x + e)*a*c^3 + 72*(2*f*x + 2*e - sin(2*f*x + 2*e))*b*c^2*d + 72*(2*f*x + 2*e - sin(2*f*x + 2*e))*a
*c*d^2 + 96*(cos(f*x + e)^3 - 3*cos(f*x + e))*b*c*d^2 + 32*(cos(f*x + e)^3 - 3*cos(f*x + e))*a*d^3 + 3*(12*f*x
+ 12*e + sin(4*f*x + 4*e) - 8*sin(2*f*x + 2*e))*b*d^3 - 96*b*c^3*cos(f*x + e) - 288*a*c^2*d*cos(f*x + e))/f

________________________________________________________________________________________

Fricas [A]
time = 0.37, size = 150, normalized size = 0.88 \begin {gather*} \frac {8 \, {\left (3 \, b c d^{2} + a d^{3}\right )} \cos \left (f x + e\right )^{3} + 3 \, {\left (8 \, a c^{3} + 12 \, b c^{2} d + 12 \, a c d^{2} + 3 \, b d^{3}\right )} f x - 24 \, {\left (b c^{3} + 3 \, a c^{2} d + 3 \, b c d^{2} + a d^{3}\right )} \cos \left (f x + e\right ) + 3 \, {\left (2 \, b d^{3} \cos \left (f x + e\right )^{3} - {\left (12 \, b c^{2} d + 12 \, a c d^{2} + 5 \, b d^{3}\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{24 \, f} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*sin(f*x+e))*(c+d*sin(f*x+e))^3,x, algorithm="fricas")

[Out]

1/24*(8*(3*b*c*d^2 + a*d^3)*cos(f*x + e)^3 + 3*(8*a*c^3 + 12*b*c^2*d + 12*a*c*d^2 + 3*b*d^3)*f*x - 24*(b*c^3 +
3*a*c^2*d + 3*b*c*d^2 + a*d^3)*cos(f*x + e) + 3*(2*b*d^3*cos(f*x + e)^3 - (12*b*c^2*d + 12*a*c*d^2 + 5*b*d^3)
*cos(f*x + e))*sin(f*x + e))/f

________________________________________________________________________________________

Sympy [B] Leaf count of result is larger than twice the leaf count of optimal. 386 vs. \(2 (165) = 330\).
time = 0.26, size = 386, normalized size = 2.26 \begin {gather*} \begin {cases} a c^{3} x - \frac {3 a c^{2} d \cos {\left (e + f x \right )}}{f} + \frac {3 a c d^{2} x \sin ^{2}{\left (e + f x \right )}}{2} + \frac {3 a c d^{2} x \cos ^{2}{\left (e + f x \right )}}{2} - \frac {3 a c d^{2} \sin {\left (e + f x \right )} \cos {\left (e + f x \right )}}{2 f} - \frac {a d^{3} \sin ^{2}{\left (e + f x \right )} \cos {\left (e + f x \right )}}{f} - \frac {2 a d^{3} \cos ^{3}{\left (e + f x \right )}}{3 f} - \frac {b c^{3} \cos {\left (e + f x \right )}}{f} + \frac {3 b c^{2} d x \sin ^{2}{\left (e + f x \right )}}{2} + \frac {3 b c^{2} d x \cos ^{2}{\left (e + f x \right )}}{2} - \frac {3 b c^{2} d \sin {\left (e + f x \right )} \cos {\left (e + f x \right )}}{2 f} - \frac {3 b c d^{2} \sin ^{2}{\left (e + f x \right )} \cos {\left (e + f x \right )}}{f} - \frac {2 b c d^{2} \cos ^{3}{\left (e + f x \right )}}{f} + \frac {3 b d^{3} x \sin ^{4}{\left (e + f x \right )}}{8} + \frac {3 b d^{3} x \sin ^{2}{\left (e + f x \right )} \cos ^{2}{\left (e + f x \right )}}{4} + \frac {3 b d^{3} x \cos ^{4}{\left (e + f x \right )}}{8} - \frac {5 b d^{3} \sin ^{3}{\left (e + f x \right )} \cos {\left (e + f x \right )}}{8 f} - \frac {3 b d^{3} \sin {\left (e + f x \right )} \cos ^{3}{\left (e + f x \right )}}{8 f} & \text {for}\: f \neq 0 \\x \left (a + b \sin {\left (e \right )}\right ) \left (c + d \sin {\left (e \right )}\right )^{3} & \text {otherwise} \end {cases} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*sin(f*x+e))*(c+d*sin(f*x+e))**3,x)

[Out]

Piecewise((a*c**3*x - 3*a*c**2*d*cos(e + f*x)/f + 3*a*c*d**2*x*sin(e + f*x)**2/2 + 3*a*c*d**2*x*cos(e + f*x)**
2/2 - 3*a*c*d**2*sin(e + f*x)*cos(e + f*x)/(2*f) - a*d**3*sin(e + f*x)**2*cos(e + f*x)/f - 2*a*d**3*cos(e + f*
x)**3/(3*f) - b*c**3*cos(e + f*x)/f + 3*b*c**2*d*x*sin(e + f*x)**2/2 + 3*b*c**2*d*x*cos(e + f*x)**2/2 - 3*b*c*
*2*d*sin(e + f*x)*cos(e + f*x)/(2*f) - 3*b*c*d**2*sin(e + f*x)**2*cos(e + f*x)/f - 2*b*c*d**2*cos(e + f*x)**3/
f + 3*b*d**3*x*sin(e + f*x)**4/8 + 3*b*d**3*x*sin(e + f*x)**2*cos(e + f*x)**2/4 + 3*b*d**3*x*cos(e + f*x)**4/8
- 5*b*d**3*sin(e + f*x)**3*cos(e + f*x)/(8*f) - 3*b*d**3*sin(e + f*x)*cos(e + f*x)**3/(8*f), Ne(f, 0)), (x*(a
+ b*sin(e))*(c + d*sin(e))**3, True))

________________________________________________________________________________________

Giac [A]
time = 0.48, size = 152, normalized size = 0.89 \begin {gather*} \frac {b d^{3} \sin \left (4 \, f x + 4 \, e\right )}{32 \, f} + \frac {1}{8} \, {\left (8 \, a c^{3} + 12 \, b c^{2} d + 12 \, a c d^{2} + 3 \, b d^{3}\right )} x + \frac {{\left (3 \, b c d^{2} + a d^{3}\right )} \cos \left (3 \, f x + 3 \, e\right )}{12 \, f} - \frac {{\left (4 \, b c^{3} + 12 \, a c^{2} d + 9 \, b c d^{2} + 3 \, a d^{3}\right )} \cos \left (f x + e\right )}{4 \, f} - \frac {{\left (3 \, b c^{2} d + 3 \, a c d^{2} + b d^{3}\right )} \sin \left (2 \, f x + 2 \, e\right )}{4 \, f} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*sin(f*x+e))*(c+d*sin(f*x+e))^3,x, algorithm="giac")

[Out]

1/32*b*d^3*sin(4*f*x + 4*e)/f + 1/8*(8*a*c^3 + 12*b*c^2*d + 12*a*c*d^2 + 3*b*d^3)*x + 1/12*(3*b*c*d^2 + a*d^3)
*cos(3*f*x + 3*e)/f - 1/4*(4*b*c^3 + 12*a*c^2*d + 9*b*c*d^2 + 3*a*d^3)*cos(f*x + e)/f - 1/4*(3*b*c^2*d + 3*a*c
*d^2 + b*d^3)*sin(2*f*x + 2*e)/f

________________________________________________________________________________________

Mupad [B]
time = 8.08, size = 183, normalized size = 1.07 \begin {gather*} \frac {2\,a\,d^3\,\cos \left (3\,e+3\,f\,x\right )-6\,b\,d^3\,\sin \left (2\,e+2\,f\,x\right )+\frac {3\,b\,d^3\,\sin \left (4\,e+4\,f\,x\right )}{4}-18\,a\,d^3\,\cos \left (e+f\,x\right )-24\,b\,c^3\,\cos \left (e+f\,x\right )-72\,a\,c^2\,d\,\cos \left (e+f\,x\right )-54\,b\,c\,d^2\,\cos \left (e+f\,x\right )+24\,a\,c^3\,f\,x+9\,b\,d^3\,f\,x+6\,b\,c\,d^2\,\cos \left (3\,e+3\,f\,x\right )-18\,a\,c\,d^2\,\sin \left (2\,e+2\,f\,x\right )-18\,b\,c^2\,d\,\sin \left (2\,e+2\,f\,x\right )+36\,a\,c\,d^2\,f\,x+36\,b\,c^2\,d\,f\,x}{24\,f} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a + b*sin(e + f*x))*(c + d*sin(e + f*x))^3,x)

[Out]

(2*a*d^3*cos(3*e + 3*f*x) - 6*b*d^3*sin(2*e + 2*f*x) + (3*b*d^3*sin(4*e + 4*f*x))/4 - 18*a*d^3*cos(e + f*x) -
24*b*c^3*cos(e + f*x) - 72*a*c^2*d*cos(e + f*x) - 54*b*c*d^2*cos(e + f*x) + 24*a*c^3*f*x + 9*b*d^3*f*x + 6*b*c
*d^2*cos(3*e + 3*f*x) - 18*a*c*d^2*sin(2*e + 2*f*x) - 18*b*c^2*d*sin(2*e + 2*f*x) + 36*a*c*d^2*f*x + 36*b*c^2*
d*f*x)/(24*f)

________________________________________________________________________________________